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ABSTRACT 


r ./The synthesis of autopilots for rotary-wing VTOL (vertical takeoff 
and landing) aircraft is more complicated than for conventional air- 
craft due to the added degrees of freedom of the rotor. In previous 
reports (see References) we have shown that improved aircraft response 
to disturbances can be obtained by considering rotor dynamics when de- 
signing the autopilot. These conclusions were based on feeding back 
perfect information on all rotor and fuselage state variables (angles, 
angular velocities and translation velocities). This report considers 
the practical situation where imperfect information on only a few state 
variables is available. Filters are designed to estimate all the state 
variables from noisy measurements of fuselage pitch/roll angles and 
from noisy measurements of both fuselage and rotor pitch/roll angles. 
The mean square response of the vehicle to a very gusty, random wind is 
computed using various filter/controllers and is found to be quite 
satisfactory although, of course, not so good as when one has perfect 
information (idealized case). 

The second part of the report considers precision hover over a 
point on the ground. A vehicle model without rotor dynamics is used 
and feedback signals in position and integral of position error are 
added. The mean square response of the vehicle to a very gusty, random 
wind is computed, assuming perfect information feedback, and is found 
to be excellent. The integral error feedback gives zero position error 
for a steady wind, and smaller position error for a random wind. 
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NOMENCLATURE 


C 

F 

G 

H 

K 

P 

Q 

R 

X 


P 

q 

u,v 


u 


x,y 


matrix of state weighting coefficients (i.e., weighting 
of e p = A e > 


matrix of control gains 

open loop (controls fixed) dynamics matrix 

control distribution matrix 

measurement scaling matrix 

matrix of Kalman filter gains 

covariance of error estimate 

power spectral density of process noise 

power spectral density of measurement noise 

mean square response matrix of state responses 

mean square response matrix of state estimates 

roll rate 

pitch rate 

longitudinal and lateral velocities of vehicle c.g. 
in vehicle axes 

control vector 

longitudinal and lateral position of vehicle c.g. 


T) integral of lateral position 

Q pitch angle 

lateral cyclic pitch 

longitudinal cyclic pitch 

| integral of longitudinal position 

T time constant 

c 
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roll rate 


r 


process noise distribution matrix 


Subscripts 

< >F 

( A ) 

C) 


rotor state 
fuselage state 
estimated state 
error in estimated state 


Superscripts 

( ) 1 matrix of augmented system 

T 

( ) matrix transpose (rows and columns interchanged) 
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1. SYNTHESIS OF AUTOPILOTS TO MAINTAIN NEARLY ZERO VELOCITY IN HOVER 
USING FUSELAGE ANGLE MEASUREMENTS WITH AND WITHOUT ROTOR ANGLE 
MEASUREMENT 

1.1 Introduction 

In a previous report [BR1] , we discussed the synthesis of hover 
autopilots by a new, more efficient, quadratic synthesis technique. We 
assumed, in that work, that perfect measurements of ali the vehicle 
states were available. In a subsequent report [BR2] we gave a detailed 
description of this new technique for "Optimal Control and Filter Syn- 
thesis by Eigenvector Decomposition. " In this report, we use the eigen- 
vector decomposition technique to synthesize filters to estimate the 
vehicle states for the tenth order system of [BR1] from measurements 
containing random errors. These filters are then combined with the 
"perfect-state-information controller" developed in [BR1] to form a 
dynamic compensating network. Finally, the mean-square response of the 
controlled vehicle to a random horizontal wind is determined. 

1. 2 A Tenth Order Model of Roll-Pitch-Horizorital Velocity for a Rotary- 

Wing VTOL Aircraft Near Hover 

A two-rigid-body model (TRBM) of a typical rotary-wing VTOL air- 
craft near hover was derived in [HA1]. The fuselage is regarded as one 
rigid body and the spinning rotor is modeled as another (axially symme- 
tric) rigid body that can be tilted with respect to the fuselage by the 
controls (longitudinal and lateral cyclic pitch) (see Fig. 1). This 
leads to a 16th order mathematical model since there are 3 translational 
degrees of freedom (DOF) for the mass center, 3 angular DOF for the 
fuselage, and 2 angular DOF for the rotor (assuming constant rotor spin 
velocity). If we disregard horizontal position of the mass center, the 
model reduces to 14th order; near hover, vertical motion and yaw motion 
are both virtually uncoupled, and the model reduces to one 10th order 
system and two uncoupled 2nd order systems (collective pitch controlling 
the vertical motion and tail rotor controlling the yaw motion). 
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FUSELAGE- FIXED AXIS FUSELAGE-FIXED 



Fig. 1-1. MODEL REFERENCE AXES. 

Our attention in this section, as in [BR1], is focused on the 10th 
order system. Using the results of many previous investigators for 
rotor- fuselage aerodynamic forces and torques, we constructed the fol- 
lowing constant-coefficient linear model for the roll-pitch-horizontal 
translation motions (linearized about equilibrium hover): 

x = Fx + Gu (1) 

where 

x T = (0 R ' < ? R ' q R 'P R ' e F ’ c P F > q F 'P F > u ' v) = state vector ’ 

T 

u = (9 ,0 ) = control vector, 

— c s 

and 
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(u,v) are velocity components of the mass center along the (x,y) fuse- 
lage axes. (0 ,0 ) are longitudinal and lateral cyclic pitch perturba- 
s c 

tions, respectively. A feedback controller, u = Cx , was designed in 
[BR1] for this system, assuming perfect measurements of the ten vehicle 
state variables. 

In Table 1, the elements of the 10x10 matrix F , the 10x2 matrix 

G , and the 2x10 matrix C , are given for the Sikorsky S-61 helicopter. 

For controls fixed (0 =0 = 0) , the characteristic roots of this 

s c 

model are quite close to those of a model used by Sikorsky engineers 
(see [BR1] ) . 


1.3 A Wind Model 


The wind components (u^.v^) along the (x,y) fuselage axes are 

modeled by independent exponentially-correlated Gauss-Markov processes 

(cf. Mil). We took the components of the wind (u and v ) to be 

w w 

independent Gauss-Markov processes defined by: 


u 

w 



+ q 


u 


1 

v = v + q 

w T w v 

c 


( 2 ) 


where 


E[q u ] = E[q v ] = 0 


and 


q (t)q (t ’ ) , q (t)q (t') _ 


2 

o 

_J 

u u u v 


w 

q v < t)q u (t,) ’ q v (t)q v (t,) _ 


-° V 


2a 


Q = 

w 


w 


x is the correlation time of the wind. In the statistical steady 


state Eq. (2) gives E[u ] = E[v ] = 0 , 

w w 

= 0 . For the simulations here, we took 

-1 


E[u 2 ] = E[v 2 ] = a 2 , E[u v ] 
w w w w 


w 


x =3.2 sec 
c 


and 


a = 20 ft 
w 


sec , which is a very strong, gusty wind. 
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Table 1-1 - Elements of the Open- Loop Dynamics Matrix, the Control 
Distribution Matrix, and a Feedback Gain Matrix for a 
Mathematical Model of the Sikorsky S-61 Helicopter 

x = Fx + Gu ; x T = t 0 R »'P R » < 3 R »P R » e F > c P F > c l F >P F > u » v ^ 

T 

u=Cx : u = [0 ,6 ] = lateral, longitudinal cyclic pitch 

— — — c s 

9 = pitch angle, cp = roll angle, ( ) = fuselage, ( ) - rotor 

r K 


u,v are longitudinal, lateral velocity components 
All units in feet, seconds; angles in radians. 


0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

- 41.3 

-601 

- 30.2 

- 42.6 

0 

0 

- 30.4 

- 50.1 

. 126 

-.284 

599 

- 56.7 

42.6 

- 29.4 

0 

0 

50.2 

- 30.2 

- . 283 

-.122 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

4 . 97 

-.94 

-.044 

.0034 

0 

0 

-.0521 

.0281 

. 00124 

-.00020 

3.53 

18.7 

-.013 

-.166 

0 

0 

-.105 

-.196 

-.00076 

-.00467 

- 15.0 

21.9 

1.03 

-.045 

- 32 . 

2 0 

4.37 

1.44 

-.0166 

.007 2 

21.9 

15.0 

-.045 

- 1.03 

0 

32 . 2 

1.44 

- 4.37 

-.0072 

-.0166 

r ° 

0 

-601 

5.52 

0 

0 

-.938 

- 4.97 

21.8 

- 16.8 “| 

L o 

0 

- 1.47 

-599 

0 

0 

1.32 

- 3.52 

- 16.8 

- 21.8 J 

r .26 

, -31 , 

.007 , 

-. 001 , 

.27 

, .97 

, .112 ,. 

. 272 , . 

000083 , 

-.0000581 

[ j .16 

, .17 , 

. 001 , 

.004 , 

-.97 

, -27 

,-. 535 , 

, 041 ,-. 

000031 , 

-.00007 9j 


These gains correspond to the performance index, 


J = 4 /o (0 F 4<P F'^! +e c )4t 
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Equation (2) is called a "shaping filter" for the wind; using this 

approach, u and v become additional (eleventh and twelfth) state 
w w 

variables of the system. They enter into Eqs. (l) only through the 
aerodynamic terms where u is replaced by u + U w and v is replaced 
by v + v^ . The augmented equations are then: 


— — 


X 


u 

w 

II 

V 

w 



u 

“ 1/T c 

0 


0 

-1/t 


v — ' 





— — 


— 


— — 


X 


G 


0 


u 

+ 

0 

u + 

1 0 


w 



— 



V 


0 


0 1 


w 


— 


_ _ 


3 W 


J V — / 
x’ 


' — ' 
G’ 


(3) 


where f u = F. >g , f y = F. 1Q , (i = 1,...10) . 

The perfect-information-controller now contains four additional 

feedback gains on u and v : 

w w 


"e ~ 


u 

c 

= Cx + C 

w 

0 

— w 

V 

L-.S-J 


L. w-J 


(4) 


The other twenty gains, C , are not changed by addition of these two 

state variables, since Eq. (2) is coupled only one-way to Eq. (1). In 

fact, the four additional gains, C , may be determined in terms of 

w 

t and the steady state Riccati matrix corresponding to Eq. (1), which 
c 

is shown in Appendix A. For the numerical values used here, the four 
additional gains are: 


C 

w 


.00021 , -.00051 
-.00050 , -.00020 


(5) 


1.4 Measurement Error Models 

Measurements of fuselage roll angle, cp , and fuselage pitch 

r 

angle, 0 , have been used for some time as basic inputs to autopilots, 

r 

They are readily obtained from the vertical reference system of the air- 
craft which usually consists of a two-degree-of-f reedom gyro with two 
electrolytic bubble levels. We have assumed that these measurements 
contain additive white noise with zero mean. For one system we assumed 
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a power spectral density of 2. 8 X 10 rad sec , which was estimated by 
considering an RMS error of .22° with a correlation time of 0.1 second. 
For another, more accurate, system, we assumed an IMS error of .09°, 
still with a correlation time of 0.1 second. Actually, vertical gyro 
measurements are sensitive to lateral accelerations since the bubble 
levels measure the "apparent vertical. " Thus our white noise error 
model is really more appropriate to angle measurements from an inertial 
measurement unit (IMU) with a stable platform. 

Measurements of the roll angle, cp , and the pitch angle, 0 , of 

the rotor tip-path-plane are more difficult and more expensive. Direct 
rotor angle measurements for articulated rotors have been made with 
potentiometers or strain gauges on the blade flapping hinge. The prin- 
cipal drawback of this technique is the required resolution of the 
measurements from rotating to non-rotating axes. This complexity and the 
extra weight are justified only for flight tests and not for operational 
use. Other direct methods could be based on directly linking the blade 
tips to a shaft mounted sensor or on electromagnetic (possibly optical) 
techniques. For example, small electromagnetic radiators at the rotor 
tips could be detected at the hub by direction-sensitive receivers, two 
in pitch and two in roll. 

Indirect rotor angle measurements have been made by sensing hub 
moments. The AH-56 helicopter originally used a shaft mounted gyro 
which was torqued by moments proportional to the hingeless blade flap- 
ping. The difficulties of this system lie in the added complexity, the 
poor drag characteristics of the gyro assembly, and in the problem of 
isolating the moment sensors from vibrations. 

We have assumed that measurements of cp and 0 contain additive 

K R 

white noise with zero mean, and we have considered two levels of power 

-6 2 

spectral density. One of these levels is high (7.1 x 10 rad sec), 

corresponding to inaccurate measurements, and the other is low (7.1 x 
—8 2 

10 rad sec), corresponding to very accurate measurements. These levels 
were estimated by considering a correlation time of 0. 1 second and RMS 
errors of 1.0 deg and 0.1 deg, respectively. The latter is probably 
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more accurate than can be expected but serves to tell us whether or not 
an attempt to develop very accurate measurements would be worthwhile in 
terms of improving vehicle response. 

We did not consider the use of measurements of the fuselage angular 

velocities, p and q . They are easily obtained with rate gyros but 
F F 

would have to be supplemented either with angle measurements (<$>,0 ) 

r r 

or velocity measurements. Measurements of fuselage velocity (u,v) 
could be obtained directly from a Doppler system or indirectly from an 
IMU. Airspeed measurements are not very accurate at the low velocities 
corresponding to hover and, of course, are not measurements of ground 
speed but speed with respect to the air. The combination of velocity 
measurements from a Doppler system with rate gyro measurements would be 
an interesting one to consider. 


1.5 A Filter to Estimate Vehicle States From Measurements of Fuselage 
Roll/Pitch Angles 

Two steady-state Kalman filters corresponding to the twelfth order 
model of Eq. (3) were designed using the eigenvector decomposition tech- 
nique of [BR2] for measurements of 0 and cp , assumed to contain 
additive white noise. Numerical values for F’ were obtained from 

Table 1 and below Eq. (2); two values for the power spectral density of 

-6 2 

the measurement noise were used, 2.8 X 10 rad sec for Filter A, and 
-6 2 

.48 x 10 rad sec for Filter B. The filters have the form 


x' = F'x' + G'u + K(z-Hx’) 


( 6 ) 


where 


z = 


+ noise 
F 

cp F + noise 


= measurements, H = 


000010000000 

000001000000 


<$’) T 


[ e R ,C PR' q R> Pp’ u ' v , u„, v w ] = estimates of the 

RRRRFFTF w.w state variables. 


7 




The gains, K , are given in Table 2, the eigenvalues of the filter 
estimates, x' = (x’-x 1 ) are given in Table 3, and the RMS errors in 
the filter estimates are given in Table 4. 


Some of the gains change sign in going from Filter A to Filter B. 
Such sign changes with measurement accuracy are not uncommon. 


The eigenvalues of the filter error, Table 3, indicate the response 
decay of the error relative to the uncontrolled and controlled modes of 
the vehicle. The smallest eigenvalues are real and indicate rather slow 
error decay for that filter mode. 

The RMS values of the estimate errors are given in Table 4. These 
are the square roots of the diagonal terms of the error covariance ma- 
trix. The effect of these errors on the vehicle response is discussed 
in Section 1.7. The rotor estimate errors and the wind velocity errors 
are not significantly changed by increased fuselage measurement accuracy, 
but the fuselage estimate errors are substantially reduced. 



Table 1-2 ,- Gain Matrices for Filters A and B 
(All units in seconds, feet; 
angles in radians) 


Filter A - Less Accurate Measurements 
.12 -3.02 -.80 5.75 .10 16.90 .22 -86.1 5.35 

1.47 1.80-6.54 .108.44 1.1535.60 2.8937.10 


7550. 1320. 

1430. -8060. 


Filter B - More Accurate Measurements 


7.60 .71 -7.40 -2.08 8.0 .20 32.2 2.0 -172. 3.56 19450. 
.19 5.36 3.30 -17.8 .20 11.92 1.91 71.3 -2.82 89.4 3510. 


3360. 

-20500 




1. 6 A Filter to Estimate Vehicle States from Measurements of Fuselage 
Roll/Pitch Angles and Rotor Roll/Pitch Angles 

Two other steady-state Kalman filters were designed for the model 

of Eq. (3); these filters have four measurements as inputs, 0 D »V 0 r> 

R R F 

and cp . Both filters assumed the less accurate fuselage angle measure- 

F -6 2 
ments (noise power spectral density 2. 8 X 10” rad sec). Two values for 

the noise power spectral density of the rotor angle measurements were 

-62 -82 
used, 7.1 x 10 rad sec for Filter C, and 7.1 x 10 rad sec for 

Filter D. The filters have the same form as (6) except that now 
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The gains K are given in Table 5, the eigenvalues of the error 
estimates, x' = x* - x' , are given in Table 6, and the RMS errors of 
the estimates are given in Table 7. 

A comparison of the RMS errors of the estimates for Filter A 
(Table 4) with those of Filter C (Table 7) shows only a slight reduc- 
tion by adding the less accurate rotor measurements. 

When the more accurate rotor measurements are used (Filter D), sig- 
nificant changes in the eigenvalues of the estimate errors occur (Table 
6). The gains on the rotor states for this Filter D are significantly 
increased (Table 5). Reduction of the rotor state estimate errors is 
accompanied by reduction of the fuselage angle and wind velocity estimate 
errors, but the velocity errors are not reduced at all (Table 7). 


Table 1-5 - Gains for Filters C and D 
(All units in seconds, 
feet; angles in radians) 


Filter C - Less Accurate Rotor Measurements 


1.99 

.13 

2. 28 

CO 

00 

to 

.47 

.08 

3.68 

.17 

-5. 10 

.86 

4560. 

610. 

.13 

1.21 

-3.42 

.96 

.002 

.42 

.28 

5.90 

-.07 

-2.32 

635. 

-3160. 

1.16 

.006 

-1.62 

-.50 

4.51 

.01 

10.5 

-.26 

-66.5 

3.38 

3890. 

704. 

.19 

1.06 

1.67 

-5.26 

.01 

7.71 

.54 

30.2 

3.91 

36.2 

794. 

-6650. 


Filter D - More Accurate Rotor Measurements 


20.9 

.05 

218. 

154. 

.27 

.40 

10.2 

2.68 

119.7 

-38.2 

57000. 

-9510. 

.05 

19.8 

-152. 

196. 

-.10 

.53 

-.67 

34.9 

-36.8 

-124. 

-9510. 

-56500. 

.007 

-.003 

-.10 

-.04 

1.68 

-.006 

1.42 

-.12 

00 

CO 

CO 

1 

.92 

392. 

4.32 

.001 

.013 

.04 

-.56 

-.006 

3.63 

.09 

6.62 

1.73 

29. 2 

-75.6 

-1117. 






Table 1-6 - Eigenvalues (sec - ^) of Estimate Errors for Filters C and D 


-13.2+5.6 j 

-1.7S±2.42j 

-4. 5S± . 31 j 

-.27+4. 10 j 

-. 0011 '!' 

-22.2±1.21j 

-12. 1+10. 9j 

-1.49+1. 42j 

-.76+ . 74j 

-.OOll 1- 


r Double roots. 


Table 1-7 - RMS Errors of Estimates for Filters C and D 
(Units in feet, seconds; angles in degrees) 





































1.7 Mean Square Response of Controlled Vehicle to a Random Wind 


Combining any one of the filters of Tables 2 or 5 with the feedback 
gains of the perfect information controller of Table 1, produces a 
dynamic compensator, shown schematically below: 

Measurement Dynamic 



Fig. 1-2. BLOCK DIAGRAM OF FIL- x = augmented state vector 
TER-CONTROLLER-AS A z = measurement signal 

DYNAMIC COMPENSATOR u = control signal 

T 

The steady-state, mean-square response of the vehicle, X = E(xx ) , 

using such an autopilot can be predicted by solving a set of 66 = ^(11 )• 

a /\y\ r T 

(12) linear equations for the elements of X = E(xx ) and adding this 
matrix to P = E(x-x)(x-x) T , * i.e., 


X = X + P 


(7) 


where 

(F-GC)X + X(F-GC) T = -KRK T (8) 

and F is the augmented (11x11) open-loop dynamics matrix, G is the 

augmented (11x2) control distribution matrix, C is the augmented 

(2x11) feedback gain matrix, K is the augmented filter gain matrix 

(11x2 or 11x4) , R is the measurement error power spectral density 

T 

matrix = E(z-Hx) (z-Hx) , and P is the covariance matrix of the 
estimate error. In addition, the steady-state mean-square control 
activity can be predicted from 

E(uu T ) = CXC T . (9) 

These computations were done using the OPTSYS computer program 
described in [BR2] and the root-mean-square (RMS) responses of the ten 


*See, e.g, , [BR3 ] , pp. 416-418. 





vehicle states and the two controls* are shown in Table 8 for autopilots 
using Filters A, B, C, and D. For comparison the RMS responses of the 
perfect information controller are also shown in Table 8; in this latter 
case, we assume perfect measurements of all twelve state variables; and 
the mean square response is obtained from 

(F-GC)X + X(F-GC) = -rQ T T (11) 

w 

E(uu T ) = CXC T (12) 

where F is defined in Eq. (3), and = power spectral density of 

q^ , defined in Eq. (2). 

As expected, the RMS responses are smallest when perfect information 
on all states is fed back; however, this is an idealized situation we 
can never realize. Using Filter A in the autopilot (the less accurate 
measurements of fuselage roll/pitch angles) still produces surprisingly 
small responses to such a strong, gusty lateral wind. 

Using Filter B in the autopilot (the more accurate measurements of 
fuselage roll/pitch angles) shows some reduction in RMS responses over 
the use of Filter A, as, expected. The interesting point is that this 
system is better than the system using Filter C which uses the less ac- 
curate measurements of both fuselage and rotor roll/pitch angles. In 
other words, Filter B estimates the rotor states more accurately than 
they are measured for Filter Cl 

Filter D uses limiting, almost unattainably accurate measurements of 
rotor pitch and roll angles along with the less accurate measurements of 
fuselage pitch and roll angles. Using Filter D in the autopilot pro- 
duces a substantial improvement in RMS response over the use of Filter A. 


* These are the square roots of the diagonal elements of 


T 

X and E (uu ) . 
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Table 1-8 - Predicted RMS Responses of Vehicle to 

Random. Wind Using Various Filter/Controllers 
(Units in feet, seconds; angles in degrees) 


S R 

^R 

q R 

P R 

S F 

*F 

q F 

P F 

u 

V 

e 

c 

e 

s 

Filter A 

.69 

.59 

7.21 

7.62 

.71 

.73 

1.28 

2.30 

3.15 

2.69 

.89 

. 95 

Filter B 

.54 

.47 

5.79 

6. 10 

.47 

.48 

.85 

1.71 

2.35 

2.01 

.80 

.84 

Filter C 

.53 

.53 

6.06 

6.78 

.55 

.66 

.85 

2.13 

2.67 

2.56 

.85 

.84 

Filter D 

. 19 

. 22 

3.19 

3.49 

.20 

. 29 

.21 

.85 

1.91 

1.96 

.68 

.65 

Perfect 

Info. 

.08 

.08 

. 16 

. 19 

.05 

.05 

.04 

.04 

.84 

.82 

.65 

.65 

Filter A - Less accurate fuselage measurements. 

Filter B - More accurate fuselage measurements. 

Filter C - Less accurate fuselage measurements and less accu- 
rate rotor measurements. 

Filter D - Less accurate fuselage measurements and more accu- 
rate rotor measurements. 

Perfect Info. - Perfect measurements of all states (limiting, 
unrealistic case). 

RMS longitudinal and lateral winds 20 ft. sec. -1 with 3.2 sec. 
correlation time. 
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2. SYNTHESIS OF AUTOPILOTS TO HOVER OVER A POINT ON THE GROUND 
2. 1 Introduction 

In Part 1, the object of the controller was to keep the vehicle 
erect (roll and pitch angles small) and keep the rotor tilt angles small. 
In this part of the report, we design controllers that go further and 
keep the vehicle close to a specific point on the ground (longitudinal 
and lateral displacements small) in the presence of a random wind. A 
well known controller design technique for constant (but unknown) input 
disturbances is "integral control"; we combine this technique here with 
quadratic synthesis, and show that it is also useful when the input dis- 
turbances (wind in this case) are exponentially-correlated random pro- 
cesses. 


2. 2 A Sixth Order Model of Roll-Pitch-Horizontal Velocity for a Rotary- 

Wing VTOL Aircraft Near Hover 

The controller designs here are based on a sixth order roll-pitch- 
horizontal velocity model of [BR1], which treats rotor-tilting as in- 
stantaneous. Table 1 gives the 6x6 open-loop dynamics matrix, the 
6x2 control-distribution matrix, and the open-loop eigenvalues of the 
model. 

This model will be augmented here by four more states x,y,|,r) de- 
fined by 




x = u , 

i - x , 



y = v , 

-=3* 

II 

where 

(x,y) 

are longitudinal and 

lateral 

center 

from a 

specified hover point 

on the 


( 1 ) 


grals of these displacements. 


As in Part 1, first-order shaping filters are used to simulate in- 
dependent random wind components, which raises the system order to 8. 
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Table 2-1 - Roll-Pitch-Horizontal Velocity Model of Sikorsky 
S-61 Helicopter for Instantaneous Rotor Tilting 


0 = fuselage pitch angle, cp = fuselage roll angle; (u,v) = 

(longitudinal, lateral) velocity components; = ^ a_ 

teral, longitudinal) cyclic pitch = (-cp ,8 ) = (-roll, pitch) 

R R 

tilt angles of rotor. All units in feet, seconds; angles in 
radians. 


— — 


— 





— 


— n 


— 

— 

0 F 


0 

0 

1 

0 

0 

0 


S F 


0 

0 

^F 


0 

0 

0 

1 

0 

0 


^F 


0 

0 

q F 


0 

0 

-.415 

.318 

.00338 

.00116 


q F 


-. 295 

6. 27 

P F 

= 

0 

0 

-1. 23 

-1.58 

.00415 

-.0124 


P F 

+ 

-23.1 

-1.08 

u 


-32.2 

0 

4.70 

-1.02 

-.0198 

-.0059 


u+u 

w 


. 977 

-32. 2 

V 


0 

32.2 

-1.02 

-4.70 

.0059 

-.0198 


v+v 
, w 


-32.2 

-.977 


Open Loop Eigenvalues are: (-1. 2, -1. 1 , . lli. 36 j , . 04±. 50j )sec 


2.3 Equilibrium State for Constant Wind; Feedback of Integral of Error 


For constant wind (u = constant, v = constant) and constant 

w w 


(0 and 9 ) , the vehicle has an equilibrium 
c s 


control deflections 

state with constant values of 9 ,cp ,u, and v , determined by the last 

four equations of Table 1 with '& = cp = 9 = cp = u = v = 0 . There 

F F F F 

are particular values of 9 and 9 that produce u = v = 0 , namely 

c s 


( 2 ) 


— — 


— — 


“ — 

e 


.00021, -.00053 


u 

c 


7 


w 

9 


-.00053 ,-.00021 


V 

s 


7 


w 


which correspond also to particular values of 0 and cp 


r — < 


- “ 


- — , 



- . 000059 . 000018 


U 

F 


7 


w 

9f _ 


.000016,. 000056 


V 

w 


(3) 


Hence, if u and v were known, additional feedback terms propor- 
w w ’ 


tional to u and v 

w w 


could, in principle, bring the vehicle to 
equilibrium with u = v = 0 . 


However, such a feedback scheme is not practical for a steady wind 


16 









since it is very sensitive to knowledge of u ,v and the elements of 

w w 

F, G, and C. Small changes (or inaccuracies) in any of these quantities 
would produce a steady-state with u and v constant but not equal to 
zero. A well-known remedy for this situation is the addition of feed- 
back proportional to the integral of the errors (in this case u and 
v). It is also quite natural in this case since the integral of velocity 
error is position, i.e., 


x = u , 

y = v . 


(4) 


[Note position (x,y) does not have to be measured to utilize this idea; 
the measured velocity components can be integrated and these "integrated 
errors" fed back. Of course, it is the measured velocity that is brought 
to zero in this case, not the true velocity.] Now, in equilibrium, x = 
y = 0=^>u = v=:0 , and the equilibrium values of the six quantities 
(0 ,cp ,x,y,6 ,0 ) are determined by the last four equations in Table 1 
and the two feedback relations: 



Thus, a "hang-off error" in (u,v) has been eliminated and there is no 

sensitivity to small changes (or inaccuracies) in u ,v , and elements 

w w 

of F, G, and C. 


Obviously, the same idea can be used again to eliminate the hang- 
off error in (x,y) if this is deemed necessary. We simply add two 
more states (|,rj) defined by 
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I = X 


t = y 


and feedback (£,ti) with four additional feedback gains C ,C ,C , ,C 

1 c£ crj s£ sr| 

In this case (x,y) have to be measured. As long as the twenty gains in 

the feedback gain matrix C are chosen so that the tenth order system is 

stable, the vehicle will come to equilibrium with x = y= u= v= 0 for 

a constant wind, with no sensitivity to small changes in u ,v ,F,G, or 

w w 

C. 


2.4 Integral Controller Design by Quadratic Synthesis 

The OPTSYS computer program was used to determine feedback gains 
for three models of the Sikorsky S-61 helicopter with and without inte- 
gral error feedback. Model A is the model of Table 1 with six states 

• • 

(6> ,cp ,0 ,cp ,u,v) . Model B is the same as A with the addition of posi- 
r r r F 

tion states (x,y) where x = u , y = v . Model C is the same as B 
with the addition of integral position states (£,t]) where | = x , 

ti = y • 


Table 2 gives the feedback gains obtained by using the following 
quadratic performance indices: 


Model A: J 


Model B: J 



F V F c 


l) 


dt 


r oo /0^-Kp 2 -t0 2 4« 2 
' o 



Model C: J 


/„ 2 2 n 2 n 2 
-co/0 -KD -t£ +0 

' r> 


+ 



dt 


where 0 =1.0 deg., x = 10.0 ft., | =50.0 ft. sec. 

o o s o 

When the longitudinal and lateral wind, u and v , are added as 

w w 

states (see Section 1.3), the gains in Table 2 are unchanged but gains 
on the wind are added. These are shown in Table 3 for the wind correla- 
tion time of 3.2 sec. 


If a filter-compensator is used with Model C, one does not "estimate" 
(|,T)) since these are simply integrals of the position estimates (or 
position measurements) (see Section 1.5). 
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2.5 Mean Square Response of Controlled Vehicle to a Random Wind 


Using the part of the OPTSYS program described in Section 1.7, 

Eqs. (11) and (12), for perfect information feedback, the RMS responses 
of the vehicle were computed for several different controllers. 

Table 4 shows the RMS vehicle responses to independent random long 
tudinal and lateral winds, each with RMS velocity of 20 ft. sec. 1 and 
correlation time 3.2 sec., for Models A, B, and C of Section 2.4, using 
the gains of Tables 2 and 3. Note the reduction in RMS position errors 
through the addition of integral control (Model B to Model C), with no 
increase in RMS control activity. 

Table 5 shows the effect of wind correlation time, t , on con- 

c 

trolled vehicle response for Model C. For « 3.2 sec., the wind is 

effectively white noise but has a very small power spectral density 
2 

(20^1^) , so it produces only small vehicle responses. For T c » 3.2 
sec., the wind is effectively constant, so with the integral control x 
and y » 0 . The largest RMS (x,y) response occurs for = 3.2 sec 

since this is in the range of time constants of. the controlled vehicle. 



























Table 2-5 - Predicted RMS Responses of Vehicle to 
Random Wind with Integral Control for 
Various Wind Correlation Times (Units in 
feet, seconds; angles in degrees) 




*F 

X 

D 

1 

B 

9 

s 

9 

c 

sk 

CO 

(steady) 

-.05 

.08 

0 

0 

-19.3 

-13.8 

-.84 

.37 

32 

.07 

.07 

.02 

.02 

.17 

.11 

.65 

.65 

3.2 

.07 

.07 

. 18 

.16 

1.3 

1.2 

.65 

.65 

. 00034 
(white) 

.02 

.02 

.04 

.04 

.25 

. 25 

.03 

.03 


Model C with perfect measurements of all states 
assumed. 

RMS longitudinal and lateral wind 20 ft. sec. - . 
*This case is calculated by solving the 10 linear 
equations x = (F-GC) -1 r w where w = v w ) T . 



















APPENDIX 


Determination of Feedback Gains on Wind 


The system equations with exponentially-decaying wind may be written 
as 

x = Fx + fw + Gu 
w = -aw 

where x = system state vector , w = wind vector , u = control vector . 
There is only a one-way coupling since x depends on w but w does 
not depend on x . 

If the performance index to be minimized is of the form 

„ 00 

IT T 

J=j\ (xAx+u Bu)dt , 

' 0 


then the Hamiltonian is 

H = £(x T Ax + u T J3u) + \ T (Fx + fw + Gu) + \ T (-aw) . 


The adjoint equations become 

. T 

X = -Ax - F X 

• T T 

\ = “ f X + a \ 

and the minimizing control is given by 

-1 T 

u = -B G \ . 

If we let 


1— * — 


r- . 


— 

X 


s , s 


X 



I w 





T | 



X 


S' a 


w 

_ wj 


L w J 


_ _ 


and differentiate this expression with respect to time, it is straight- 
forward to show that S,S w> ant * cr must satisfy the following differ- 
ential equations: 
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• T -IT 

S = -SF - F S - A + SGB G S 

S = -(F-GB -1 G T S) T S + S a - Sf 
w w w 

• T T T_ -1 T 

0 =aa+aa-fS -Sf+SGB GS 

WWW w 

The matrix Riccati equation for S is the same as if there were no 
wind. The equation for is linear but depends on S , and the equa- 

tion for a is also linear but depends on S w . The feedback law may 
be written as 

u = -Cx - C w , 
w 

where 

-IT -IT 

C = B GS, C = B GS , 
w w 

so, in principle, the gains , which depend on S^ , may be calculated 

separately after C has been determined from S . 

For the steady-state case, it is easier to use OPTSYS to compute C 
and C w together by eigenvalue decomposition. Nonetheless, it is im- 
portant to note that C is not changed by addition of a wind model. 
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